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Abstract. Solving the full classical four-body Coulomb problem numerically 
using a Wigner initial distribution we formulate a classical-quantum hybrid 
approach to study triple ionization by single photon absorption from the Li ground 
state in the threshold region. We confirm the Wannier threshold law cr oc E" and 
we show that the a determined in the interval between 2 — 5 eV deviates from the 
analytical threshold value of 2.16 which we find in the interval between 0.1 — 2 
eV. 

PACS numbers: 3.65.Sq, 32.80.Fb, 34.80.Dp 

Triple photoionization of Lithium is the most fundamental atomic process involving 
three bound electrons. Only recently Wehlitz et al. succeeded in measuring the 
triple photoionization cross section down to 2 eV above threshold [1,2]. Subsequent 
experiments have produced various double ionization cross sections, also very close to 
threshold. For Lithium [3] and most recently Beryllium [4], it was demonstrated 
convincingly that the double photoionization cross section has small oscillations 
superimposed on the rising smooth cross section. This has cast some doubt on the 
validity of Wannier's (classically derived) threshold law [5,6] which predicts very close 
to threshold a power law behavior of the cross section, 

aiE^) oc {EJI - If , (1) 

where is the photon energy, / the respective threshold energy and a a characteristic 
exponent which is related to the stability of a classical fixed point of the iV-electron 
dynamics [7]. 

For the present case of Lithium triple ionization (/ ~ 7.478 a.u.) the experimental 
results reach only down to 2 eV above threshold [2]. The corresponding fit of the 
experimental data with Eq. ^ yields acxp = 2.05 close to the theoretical value of 
a = 2.16 derived by Klar and Schlecht [8]. It is known from the well studied two- 
electron escape case [9] that a Wannier exponent fitted to a cross section over an 
energy interval which is close but a finite distance away from threshold yields always 
a smaller a than the analytically predicted one. Hence, the experimental value is 
consistent with Wannier's prediction but not conclusive since one does not know what 
happens closer to threshold. 

The present classical study deals for the first time with the full four-body problem 
close to threshold and provides the triple photoionization (TPI) probability starting at 
E — 0.9 eV excess energy. We are able to confirm that, indeed, the Wannier threshold 
law with an exponent of a = 2.16 is reached, but only for energies E = E^ — I < 2eV. 
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By successively fitting finite energy intervals above the threshold / to our result, we 
can confirm the experimental result for a. 

We formulate the TPI process from the Li ground state (ls^2s) as a two step 
process [10,11]. First, one electron absorbs the photon (photo-electron). Then, due to 
the electronic correlations, redistribution of the energy takes place resulting in three 
electrons escaping to the continuum. We express the above two step process as 

fT3+ =CTabsP'+, (2) 

where aabs is the total absorption cross section and is the probability for triple 
ionization. In what follows, we evaluate P'^+ and use the experimental data of 
Wehlitz [12] for (Tabs- Physically interpreted, this relation splits the photon absorption 
((Tabs) from the subsequent energy redistribution in the three-electron system. The 
latter can lead to TPI which we calculate in phase space formally from 

P3+ ^ lini / drp3+ exp((i - iabs)^ci)p(r), (3) 

with the classical Liouvillian given by the Poisson bracket {H, } [13] propagated 
from the time t^hs of photo absorption, with H the four-body Coulomb Hamiltonian. 
The primary electron absorbs the photon at the nucleus (ri = 0), an approximation 
which becomes exact in the limit of high photon energy [14]. We note that no account is 
taken of the direction of polarization of the incident photon, since we currently consider 
electron orbitals that are spherically symmetric. Immediately after absorption, the 
phase space distribution of the remaining two electrons is the Wigner transform of 
the corresponding initial wavefunction ip{ri = 0, r2,r3), where the are the electron 
vectors starting at the nucleus. 

In general, close to threshold the ionization probability does not depend on 
the details of the initial wavefunction [5]. Hence, we approximate it as a simple 
product of hydrogenic orbitals (j)f'{ri) with effective charges Zi to facilitate the Wigner 
transformation. The Zi are chosen to reproduce the known ionization potentials Ii, 
namely for the 2s electron Z3 = 1.259 {I3 — 0.198 a. u.) and for the Is electron 
Z2 = 2.358 {I2 = 2.780 a. u.). We use atomic units throughout the paper if not stated 
otherwise. The Wigner distribution W conserves energy only in the mean [15]. Near 
E = 0, however, energy conservation is vital. Therefore, the Wigner distributions for 
the individual electron orbitals (j^f^lji), W^z^ , are restricted to their respective energy 

shell, leading to the initial phase space distribution 

pir)^^5{ri)d{ei+h-Lo) Yl W^z,{r„p,)d{e,+I,), (4) 

1=2,3 

where ei, i — 1,2,3 are the individual electron energies. The advantages of the Wigner 
distribution as an initial phase space distribution over other distributions as well as 
the reasons for restricting the Wigner distribution of the hydrogenic orbitals on their 
respective energy shell is discussed in ref. [16]. 

With p{T) from Eq. ^ the initial phase space volume to be sampled reduces 
significantly, although regularized coordinates [17] are required to avoid problems with 
electron trajectories starting at the nucleus. Other than that, the integral in Eq. Q 
is evaluated with a standard Monte-Carlo technique which entails following classical 
trajectories in phase space (CTMC) [16,18-20]. The projector V'^^ indicates that 
we integrate only over those parts of phase space that lead to TPI. Triple ionization 
is decided by propagating trajectories long enough so that the individual electron 
energies ei, i = 1,2,3, are positive and stable within a margin that guarantees all 
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Figure 1. TPI cross section obtained by multiplying the TPI probability 
from the present calculation with the total photo cross section from [12] (o) in 
comparison to the experiment [2] (•). 



three electrons are very far away from the nucleus and each other. A similar approach 
to the one described above was successfully used to describe the knockout mechanism 
for the double ionization of He from the ground [21], the excited states [22], and 

the double ionization of H2 [23] . 

We also note that generally, not only close to threshold, the classical TPI 
trajectories provide information as to how energy is redistributed from the primary 
photo electron to the other two electrons which is analyzed in detail in [24] . 

Figure n shows the TPI cross section cr^+ resulting from the probability P^^ in 
connection with Eq. We find very good agreement with the experimental results. 
Thus, our classical approach with an approximate initial quantum wavefunction 
captures the relevant correlations among the three electrons which mainly form after 
the photo absorption at lower excess energies. Let us note that the reason we currently 
consider energies of 0.9 eV and above is the numerical difficulty involved in computing 
P^'^. Specifically, in order to obtain 10'^ TPI trajectories at E ~ 0.9 eV one has 
to evolve 10^° trajectories with the CTMC method. One may be tempted to see a 
slightly different slope of the experimental curve compared to the theoretical one in 
Figure^ Concerning the theoretical curve this may be due to the fact that we employ 
a "high energy" approximation for the photo absorption, namely that the photon is 
absorbed by an electron that is sitting initially at the nucleus (ri =0). In the future 
we may relax this approximation by using a photon-frequency dependent assumption 
for r(uj) of the photon-electron as described in ref. [25]. However, given the present 
accuracy of the experiment near threshold the slope of the experimental curve is 
somewhat uncertain as well, see error bars in Figure ^ Luckily these experimental 
and theoretical difficulties do not hamper the present goal of analyzing the behavior 
of the cross section towards threshold, E ^ 0. We want to investigate whether the 
Wannier power law for cr^^ is really approached in the limit of vanishing threshold 
energy. To this end we have fitted a^~^ = ao{E/I)°', where (Jq and a are fit parameters 
while / is the triple ionization potential. For the closest energy interval to threshold 
we could reach, 0.9 eV< E <2 eV, athoo = 2.15 very close to the analytical value of 
a = 2.16, as Figure |21 reveals. We then apply the fit, keeping the lower limit of the 
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Figure 2. The Wannier exponent a, obtained by fitting finite energy intervals 
to Eq. GJ, see text. The filled circle is the experimental value from [2]. 



energy interval constant, E = 0.9 eV, and increasing the upper energy limit until we 
reach E = 4.0 eV (acxp is obtained in the range 2 — 5.1 eV = 3.1 eV). Subsequently, 
we shift the 3.1 eV interval to higher excess energies to obtain a as a function of the 
upper limit of the energy range [26]. 

Figure El illustrates that the Wannier exponent fitted to the cross section over 
an energy interval which is close but a finite distance away from threshold yields a 
smaller a than the analytically predicted one. Hence, the fit to the experimental data 
(filled circle in Figure ^1 is consistent with Wannier's theory and in good agreement 
with our present theoretical result. The strong variation of the Wannier exponent in 
figure|21also indicates that the threshold region where Wannier's threshold law applies 
is certainly less than the energy range shown in Figure |3 

Secondly, we explore if the triple photoionization cross section for different atomic 
elements have a similar shape close to the threshold region. We use the shape formula 
for the TPI cross section [27], 



to obtain a dimensionless cross section u/ctm as a function of the dimensionless excess 
energy x — E / Em with i?M, ctm as fitting parameters and with a set to its analytical 
value of 2.16. Eq.Q reproduces, by construction, for Coulomb complete break-up 
processes, the Wannier threshold law and the cross section for high excess energies. 
In Figure 131 one sees that the experimental data for Li, Ar, and Ne fall on top of the 
theoretical data from Figure ^ 

It has been argued that a secondary power law, or at least additional structure 
in the TPI cross section could originate from the very different binding energies for 
the electrons in the Lithium atom with its ls^2s configuration in contrast to Neon 
and Argon which both contribute three electrons from a single shell to the ionization 
(2p and 3p, respectively) [2]. The present calculation does take into account the 
difference in binding energy and spatial extension of the respective orbitals. We find 
our results to be consistent with a smooth change of the TPI cross section for the 
ground state of Li. Hence, one may conclude from the agreement of our Lithium 
calculation regarding the shape of the ionization cross section with Neon and Argon 
that the different binding energies do not strongly influence the shape of the cross 
section. 

In summary in the framework of a quantum-classical hybrid approach we have 
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Figure 3. TPI cross section in scaled coordinates [27] for Lithium (this work, 
•), and from the experiments on Lithium ( [2], o), Neon ( [26], /\) and Argon 

( [26], □). 



analyzed triple photoionization of Lithium near threshold using a Wigner initial 
distribution and a classical propagation of the three electrons in time. We can 
confirm that the total cross section grows from threshold with a power of a = 2.16 in 
accordance with Wannier's threshold theory. In addition, using the shape formula, we 
find our results to be consistent with a smooth change of the triple ionization cross 
section for the Li ground state. 

We gratefully acknowledge discussions with Thomas Pattard and thank Ralf 
Wehlitz for providing his data in electronic form. 
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